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We show the construction of the dual superconducting theory for the confinement mech- 
anism from QCD in the maximally abelian (MA) gauge using the lattice QCD Monte Carlo 
simulation. We find that essence of infrared abelian dominance is naturally understood 
with the off-diagonal gluon mass m s ~ 1.2GeV induced by the MA gauge fixing. In the 
MA gauge, the off-diagonal gluon amplitude is forced to be small, and the off-diagonal 
gluon phase tends to be random. As the mathematical origin of abelian dominance for 
confinement, we demonstrate that the strong randomness of the off-diagonal gluon phase 
leads to abelian dominance for the string tension. In the MA gauge, there appears the 
macroscopic network of the monopole world-line covering the whole system. We investi- 
gate the monopole-current system in the MA gauge by analyzing the dual gluon field B^. 
We evaluate the dual gluon mass as m# = 0.4 ~ 0.5GeV in the infrared region, which is 
the lattice-QCD evidence of the dual Higgs mechanism by monopole condensation. Ow- 
ing to infrared abelian dominance and infrared monopole condensation, QCD in the MA 
gauge is describable with the dual Ginzburg-Landau theory. 

I. QCD and Dual Superconducting Theory for Confinement 

Since 1974, quantum chromodynamics (QCD) has been established as the fundamental 
theory of the strong interaction, however, it is still hard to understand the nonperturbative 
QCD (NP-QCD) phenomena such as color confinement and dynamical chiral-symmetry 
breaking, in spite of the simple form of the QCD lagrangian 

£qcd = —XjG^GT + q(i P - m q )q. (1) 

In particular, to understand the confinement mechanism is one of the most difficult prob- 
lems remaining in the particle physics. As the hadron Regge trajectory and the lattice 
QCD simulation show, the confinement force between the color-electric charges is char- 
acterized by the one-dimensional squeezing of the color-electric flux and the universal 
physical quantity of the string tension a ~ lGeV/fm. 

As for the confinement mechanism, Nambu first proposed the dual superconducting the- 
ory for quark confinement, based on the electro-magnetic duality in 1974. 1 In this theory, 
there occurs the one-dimensional squeezing of the color-electric flux between quarks by the 
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dual Meissner effect due to condensation of bosonic color-magnetic monopoles. However, 
there are two large gaps between QCD and the dual superconducting theory. 2 

(1) The dual superconducting theory is based on the abelian gauge theory subject to the 
Maxwell- type equations, where electro- magnetic duality is manifest, while QCD is a 
nonabelian gauge theory. 

(2) The dual superconducting theory requires condensation of color-magnetic monopoles 
as the key concept, while QCD does not have such a monopole as the elementary 
degrees of freedom. 

These gaps can be simultaneously fulfilled by the use of the MA gauge fixing, which reduces 
QCD to an abelian gauge theory. In the MA gauge, the off-diagonal gluon behaves as a 
charged matter field similar to in the Standard Model and provides a color-electric 
current in terms of the residual abelian gauge symmetry. As a remarkable fact in the MA 
gauge, color- magnetic monopoles appear as topological objects reflecting the nontrivial 
homotopy group n 2 (SU(A" c )/U(l) iVc ^ 1 ) =Z^ C_1 , similarly in the GUT monopole. 3 " 6 

Thus, in the MA gauge, QCD is reduced into an abelian gauge theory including both 
the electric current j M and the magnetic current k^, which is expected to provide the 
theoretical basis of the dual superconducting theory for the confinement mechanism. 

2. MA Gauge Fixing and Extraction of Relevant Mode for Confinement 

In the Euclidean QCD, the maximally abelian (MA) gauge is defined by minimizing 2,6 

R of[ [A,(.)} = J d 4 xtr[D„H][D„H^ = ^J d 4 * ]T |^(*)| 2 , (2) 

with the SU(iV c ) covariant derivative operator = + ieA^ and the Cartan decom- 
position A^(x) = A^(x) ■ H + J2 a A^(x)E a . In the MA gauge, the off-diagonal gluon 
components are forced to be as small as possible by the SU(N C ) gauge transformation. 
Since the covariant derivative obeys the adjoint gauge transformation, the local form 
of the MA gauge fixing condition is derived as 2,6 

[H,[D tl ,[D M H]]]=Q. (3) 

(For N c = 2, this condition is equivalent to the diagonalization of $ma = [D^, [D^,r 3 ]], 
and then the MA gauge is found to be a sort of the 't Hooft abelian gauge 3 .) In the MA 
gauge, the gauge symmetry G = SU(A" c ) loca i is reduced into H = U(l)^" 1 x Weyl^^, 
where the global Weyl symmetry is the subgroup of SU(iV c ) relating the permutation of 
the N c bases in the fundamental representation. 2 ' 6 

We summarize abelian dominance, monopole dominance and extraction of the relevant 
mode for NP-QCD observed in the lattice QCD in the MA gauge. 

(a) Without gauge fixing, all the gluon components equally contribute to NP-QCD, and 
it is difficult to extract relevant degrees of freedom for NP-QCD. 
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(b) In the MA gauge, QCD is reduced into an abelian gauge theory including the electric 
current j M and the magnetic current k^. The diagonal gluon behaves as the abelian 
gauge field, and the off-diagonal gluon behaves as the charged matter field. In the MA 
gauge, the lattice QCD shows abelian dominance for NP-QCD (confinement 6,7 , chiral 
symmetry breaking 8 , instantons 9 ) : only the diagonal gluon is relevant for NP-QCD, 
while off-diagonal gluons do not contribute to NP-QCD. In the lattice QCD, there 
appears the global network of the monopole world-line covering the whole system in 
the MA gauge. (See Fig.3(a).) 

(c) The diagonal gluon can be decomposed into the "photon part" and the "monopole 
part", corresponding to the separation of j M and k^. In the MA gauge, the lattice 
QCD shows monopole dominance for NP-QCD: the monopole part (k^ ^ 0, j M = 0) 
leads to NP-QCD, while the photon part (j M ^ 0, k^ = 0) seems trivial like QED and 
does not contribute to NP-QCD. 

Thus, by taking the MA gauge, the relevant collective mode for NP-QCD can be extracted 
as the color- magnetic monopole. 2 ' 10 

3. Essence of Abelian Dominance : Off-diagonal Gluon Mass in MA Gauge 

In this section, we study essence of abelian dominance for NP-QCD in the MA gauge 
in terms of the generation of the effective mass m s of the off-diagonal (charged) gluon 
by the MA gauge fixing. In the SU(2) QCD partition functional, the mass generation of 
the off-diagonal gluon A^ = (A 1 ± 2A 2 )/\/2 in the MA gauge is expressed as 2 

Z^ D = J DA,eMiS QC D[A,]}5($± A [A,])A PF [A»] 

= J DAlexp{zS cS [Al}} J £Wjexp{i J d 4 x m^A+AlL}^], (4) 

with $ MA = [Dfi, [-D M , r 3 ]], the Faddeev- Popov determinant A FP , the abelian effective 
action S^LA^] and a smooth functional 

To investigate the off-diagonal gluon mass m g, we study the Euclidean gluon prop- 
agator Gf v (x - y) = (A a ^(x)A b u (y)) in the MA gauge, using the SU(2) lattice QCD. 2 
As for the residual U(l)3 gauge symmetry, we impose the U(l)3 Landau gauge fixing 
to extract most continuous gauge configuration under the MA gauge constraint and to 
compare with the continuum theory. The continuum gluon field A a {x) is extracted from 
the link variable as U^(s) = exp(iaeA^(s)^-). Here, the scalar-type gluon propagator 
G^(r) = Y,f l =i(A^(x)A^(y)} is useful to observe the interaction range of the gluon, 

because it depends only on the four-dimensional Euclidean distance r = ^ (x^ — y^) 2 . 

We show in Fig.l(a) G^(r) and G+-(r) = K=M + ^)A~{y)) = \{G^{r) + G 2 »} 
in the MA gauge using the SU(2) lattice QCD with 2.2 < (3 < 2.4 and the various lattice 
size (12 3 x 24, 16 4 , 20 4 ). Since the massive vector-boson propagator with the mass M 
takes a Yukawa-type asymptotic form as G^(r) ~ 7375- exp(— Mr), the effective mass 
m off of the off-diagonal gluon A^(x) can be evaluated from the slope of the logarithmic 
plot of r 3 / 2 G^~ (r) ~ exp(— m fjr) as shown in Fig. 1(b). The off-diagonal gluon A^(x) 
behaves as the massive field with m g ~ 1.2 GeV in the MA gauge for r ^ 0.2 fm. 

We perform also the standard mass measurement for the off-diagonal gluon in the similar 
manner to the hadron mass measurement in the lattice QCD. We calculate 3>^(t) = 
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Figure 1. (a) The scalar-type gluon propagator G" (r) as the function of the 4- 
dimensional distance r in the MA gauge in the SU(2) lattice QCD with 2.2 < (5 < 2.4 
and various lattice size (12 3 x 24, 16 4 , 20 4 ). (b) The logarithmic plot for the scalar 
correlation r z l 2 G a „(r). The off-diagonal gluon propagator behaves as the Yukawa- 



type function, ~ <EELjjMn , ( c ) The logarithmic plot of the temporal correlation 
r+~(r) = ($+(r)$~(0)) as the function of the temporal distance r in the lattice QCD 
with 2.3 < (3 < 2.35 with 16 3 x 32 and 12 3 x 24. From the slope of the dotted lines in (b) 
and (c), the effective mass of the off-diagonal gluon can be estimated as m s ~ 1.2GeV. 



/ dxAt(x, r), and measure the temporal correlation of r+~(r) = ($+(r)$~(0)) in the MA 
gauge with the U(l) 3 Landau gauge. We obtain the off-diagonal gluon mass m oS ~ 1.2GeV 
again from the slope of the logarithmic plot of r+7(r) as the function of the temporal 
distance r in the lattice QCD with 2.3 < (3 < 2.35 with 16 3 x 32 and 12 3 x 24, as shown 
in Fig. 1(c). 

Thus, essence of infrared abelian dominance in the MA gauge can be physically inter- 
preted with the effective off- diagonal gluon mass m off induced by the MA gauge fixing. 2 
Due to the effective mass m Q g ~ 1.2GeV, the off-diagonal gluon can propagate only 
within the short range as r < m~g ~ 0.2fm, and cannot contribute to the infrared QCD 
physics in the MA gauge, which leads to abelian dominance for NP-QCD. 2 

4. Randomness of Off-diagonal Gluon Phase as the Mathematical Origin of 
Abelian Dominance for Confinement 

In the lattice QCD, the SU(2) link variable is factorized as U^(s) = M At (s)w At (s), accord- 
ing to the Cartan decomposition SU(2)/U(1) 3 x U(l) 3 . Here, u M (s) = exp{ir 3 6> 3 (s)} G 
U(l) 3 denotes the abelian link variable, and the off-diagonal factor M^(s) G SU(2)/U(1) 3 
is parameterized as 

M,(s) = e^l^em = ( cosfl (s) -am^e^M \ { ) 

MV ' \smd^(s)e lXtl(s) cos0 M (s) / v ' 

In the MA gauge, the diagonal element cos# M (s) in M^(s) is maximized by the SU(2) 
gauge transformation and the "abelian projection rate" becomes almost unity as i?Abei = 
(cos^ m (s))ma — 0.93 at (3 — 2.4. Using the lattice QCD simulation, we find the two re- 
markable features of the off-diagonal element e zx ^ s ' sin# M (s) in M M (s) in the MA gauge. 2 ' 6 
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(1) The off-diagonal gluon amplitude |sm^ M (s)| is forced to be minimized in the MA 
gauge, which allows the approximate treatment on the off-diagonal gluon phase. 

^ The off-diagonal phase variable Xm( s ) i s not constrained by the MA gauge- fixing con- 
dition at all, and tends to be random. 

Therefore, x^( s ) can be regarded as a random angle variable on the treatment of M^(s) 
in the MA gauge in a good approximation. 

Now, we show the analytical proof of abelian dominance for the string tension or the 
confinement force in the MA gauge, within the random- variable approximation for Xn( s ) 
or the off-diagonal gluon phase. 2 ' 6 Here, we use (e* XM ^^MA — fo* dx^(s) ex P{^XAt( s )} — 0. 
In calculating the Wilson loop (W C [U]) = (trU c U^s)) = (trIT c {M^(s)^(s)}), the off- 
diagonal matrix M^(s) is simply reduced as a c-number factor, M M (s) — > cos# M (s) 1, and 
then the SU(2) link variable C/ M (s) is reduced to be a diagonal matrix, 

U^s) = M^(s)u^s) -> cos^(s)^(s), (6) 

after the integration over X/i( s )- For the RxT rectangular C, the Wilson loop Wc?[C/] in 
the MA gauge is approximated as 

(W C [U\) = (trnf =1 {M^( Sl )^( Si )}) ~ (nf =1 cos^( Si )-trn^ =1 ^( S ,))MA 

~ (exp{Ef =1 ln(cos^ i (s i ))}) M A (W c [u}) MA , (7) 

with the perimeter length L = 2(R+T) and the abelian Wilson loop Wc[it] = trU^u^Si). 
Replacing Y^=\ m { cos (^U( s «))} by its average L(ln{cos0 M (s)})MA in a statistical sense, 
we derive a formula for the off-diagonal gluon contribution to the Wilson loop as 2 ' 6 

= (Wc[U])/(W c [u])ma ^ exp{L(ln{cos^(s)}) M A}, (8) 

which provides the relation between the macroscopic quantity Wff and the microscopic 
quantity (ln{cos ^(s)}}ma- Using the lattice QCD, we have checked this relation for large 
loops, where such a statistical treatment is accurate. 2,6 

In this way, the off-diagonal gluon contribution W£f obeys the perimeter law in the 
MA gauge, and then the off-diagonal gluon contribution to the string tension vanishes as 

R + T 

0"su(2) - OAbei ~ -2(m{cos6> M (s)}) M A lim — = 0. (9) 

K,l ^oo txL 

Thus, abelian dominance for the string tension, 0su(2) — 0Abeb can be demonstrated in 
the MA gauge within the random-variable approximation for the off-diagonal gluon phase. 
Also, we can predict the deviation between <7gu(2) and OAbei as °"su(2) > CAbeb d ue to the 
finite size effect on R and T in the Wilson loop. 2,6 

5. The Structure of QCD-Monopoles in terms of the Off-diagonal Gluon 

Let us compare the QCD-monopole with the point-like Dirac monopole. There is no 
point-like monopole in QED, because the QED action diverges around the monopole. 
The QCD-monopole also accompanies a large abelian action density inevitably, however, 
owing to cancellation with the off- diagonal gluon contribution, the total QCD action is 
kept finite even around the QCD-monopole. 2 ' 6 
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Figure 2. (a) The total probability distribution P(S) on the whole lattice and (b) 
the probability distribution Pk(S) around the monopole for SU(2) action density Ssu(2) 
(dashed curve), abelian action density S'Abei (solid curve) and off-diagonal gluon con- 
tribution S s (dotted curve) in the MA gauge at (3 — 2.4 on 16 4 lattice. Around the 
QCD-monopole, large cancellation between SAbei and S ff keeps the total QCD-action 
small, (c) The schematic figure for the QCD-monopole structure in the MA gauge. The 
QCD-monopole includes a large amount of off-diagonal gluons around its center. 



To see this, we investigate the structure of the QCD-monopole in the MA gauge in terms 
of the action density using the SU(2) lattice QCD. 2 From the SU(2) plaquette P^ {2) (s) 
and the abelian plaquette P^ >e[ (s), we define the "SU(2) action density" S^ 2 '(s) = 
1 - |trPy a )(s), the "abelian action density" S^ bcl (s) = 1 - |trP^ bel (s) and the "off- 
diagonal gluon contribution" S?^(s) = Sf^( 2 '(s) — S^ el (s). In the lattice formalism, the 
monopole current k^s) is defined on the dual link, and there are 6 plaquettes around the 
monopole. Then, we consider the average over the 6 plaquettes around the dual link, 

1 i 

S(s,fi) = H \^^\S a ^{s + m^f). (10) 

iZ a/3-y m=0 

We show in Fig. 2(b) the probability distribution of the action densities S SU (2), ^Abei and 
S Q g around the QCD-monopole in the MA gauge. We summarize the results on the QCD- 
monopole structure as follows. 

(1) Around the QCD-monopole, both the abelian action density S'Abei and the off-diagonal 
gluon contribution S Q g are largely fluctuated, and their cancellation keeps the total 
QCD-action density Ssu(2) small. 

(2) The QCD-monopole has an intrinsic structure relating to a large amount of off- 
diagonal gluons around its center, similar to the 't Hooft-Polyakov monopole. 

(3) At the large-distance scale, off-diagonal gluons inside the QCD-monopole become 
invisible, and the QCD-monopole can be regarded as the point-like Dirac monopole. 

(4) From the concentration of off-diagonal gluons around QCD-monopoles in the MA 
gauge, we can naturally understand the local correlation between monopoles and in- 
stantons. In fact, instantons tend to appear around the monopole world-line in the 
MA gauge, because instantons need full SU(2) gluon components for existence. 2,9 
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6. Lattice-QCD Evidence of Infrared Monopole Condensation 

In the MA gauge, there appears the global network of the monopole world-line covering 
the whole system as shown in Fig.3(a), and this monopole- current system (the monopole 
part) holds essence of NP-QCD 2 ' 8-10 . We finally study the dual Higgs mechanism by 
monopole condensation in the NP-QCD vacuum in the MA gauge. 

Since QCD is described by the "electric variable" as quarks and gluons, the "electric 
sector" of QCD has been well studied with the Wilson loop or the inter-quark potential, 
however, the "magnetic sector" of QCD is hidden and still unclear. To investigate the 
magnetic sector directly, it is useful to introduce the "dual (magnetic) variable" as the 
dual gluon field B^, which is the dual partner of the diagonal gluon and directly couples 
with the magnetic current k^. 

Due to the absence of the electric current j M in the monopole part, the dual gluon 
can be introduced as the regular field satisfying (d A B)^ u = *F fJjU and the dual Bianchi 
identity, d^*(d A B)^ = j v = 0. By taking the dual Landau gauge d^B^ = 0, the field 
equation is simplified as d 2 B ti = k^, and therefore we obtain the dual gluon field B^ from 
the monopole current k^ as 



In the monopole-condensed vacuum, the dual gluon B^ is to be massive, and hence we 
investigate the dual gluon mass m B as the evidence of the dual Higgs mechanism. 

First, we put test magnetic charges in the monopole- current system in the MA gauge, 
and measure the inter-monopole potential Vm{t) to get information about monopole con- 
densation. Since the dual Higgs mechanism is the screening effect on the magnetic flux, 
the inter-monopole potential is expected to be short-range Yukawa-type. Using the dual 
Wilson loop Wd as the loop-integral of the dual gluon, 



We show in Fig. 3(b) the inter-monopole potential V M (r) in the monopole part in the MA 
gauge. 2 Except for the short distance, the inter-monopole potential can be almost fitted 
by the Yukawa potential Vm{t) = — e ™ Br , after removing the finite-size effect of the 
dual Wilson loop. In the MA gauge, the dual gluon mass is estimated as m B ~ 0.5GeV 
from the infrared behavior of Vm{t)- 

Second, we investigate also the scalar-type dual gluon propagator (B^^B^y)) ma as 
shown in Fig.3(c), and estimate the dual gluon mass as m B ^ 0.4 GeV from its large- 
distance behavior. 

From these two tests, the dual gluon mass is evaluated as mg = 0.4 ~ 0.5 GeV, and this 
can be regarded as the lattice-QCD evidence for the dual Higgs mechanism by monopole 
condensation at the infrared scale. 




(11) 




the potential between the monopole and the anti-monopole can be derived as 
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Figure 3. The SU(2) lattice-QCD results in the MA gauge, (a) The monopole world-line 
projected into R 3 on the 16 3 x 4 lattice with (3 = 2.2 (the confinement phase). There 
appears the global network of monopole currents covering the whole system, (b) The 
inter-monopole potential Vm{t) v.s. the 3-dimensional distance r in the monopole-current 
system on the 20 4 lattice. The solid curve denotes the Yukawa potential with m B = 
0.5GeV. The dotted curve denotes the Yukawa-type potential including the monopole-size 
effect, (c) The scalar-type dual-gluon correlation ln(r E (B^xjB^y)) ma) as the function 
of the 4-dimensional Euclidean distance te on the 24 4 lattice. 



To summarize, the lattice QCD in the MA gauge exhibits infrared abelian dominance 
and infrared monopole condensation, and therefore the dual Ginzburg-Landau (DGL) 
theory can be constructed as the infrared effective theory directly based on QCD in the 
MA gauge. 2 . 
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